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An algorithm is suggested for defining the coordinants of 
X-ray and Gamma ray radiation sources by measuring the time 
lag of the transmission of radiation flares from various 
points of space i n .whtch spacecraft are located. Instances 
where the 2-x, 3-x and 4x spacecraft are used. 



The Problem of Local i za t 1 on 


I . 


A sudden and intermittent increase in the intensity of x-ray 
or gamma radiation which is detected by spacecraft is commonly re- 
ferred to as x-ray or gamma-ray flares Flares are observed 

comparatively seldom. Therefore, when registering them in n 
various points (which we will call reference points) , 1 = q 

one may subesequently indentify "temporal spectra" of^^*(i.e. the 
dependence of the radiations intensity on time), which pertain to 
one and the same flare. Using statistical methods considered in [ 2 j. 
the size of the time lagXj^ , M j of the flare's trans- 
mission in points (as compared to pointj^Q, which is taken 

for the original), may be determined by these spectra. 

The rates of time lags are used for the localization of a 
flare's source, i.e. for determining three coordinants of the 
source , . or only of the direction towards the source if it 

is presume'd that the source may be considered infinitely remote. 
Together withF.^, and ^ i t is al so. necessary to determine the statis- 
tical characterisitics of their errors ^P.and/\f ... The methods 

y ^ V 

of locating sources are similar to the methods of hyperbolic naviga- 
tion (Loran's system) [ 3 J . In view of the fact that the solution 
to the problem of localization may not turn out to be simple, it is 
advisable to use additional information on the orientation of the 
axes of the detectors which register the flares. .If the diagrams of 
the detectors' directionality are not spheri cal , then the information 
may be essential for uncovering the reason that the solution is not 
simple. 

Thus it is supposed that when locating the sources the following 
initial data is used: 

a) the rate of time lags -l^n. “~jand 

the covariance matrix of their error* 
which is determined by the spectra 
=Q 9/1 1^2]’ also the moment of Moscow 

(or Greenwich) time f of the flare's trans- 
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emission in the point M q (since the 
flare is not impulsive then corresponds 
to some representative point of the spec- 
trum So>"'ts origin, for example); 

b) the radi i -vectors /*. of the points /M .., 

{i = 0« “ I-yand the covariance matrixes 

. ^ T{ of their errors ^ at the moments 

+ r ' ( will be made to equal 
zero)' in some system of coordinants OXYZ; 

c) the single vectors of direction of the axes 

of the diagrams of the detectors' direction- 
alities di in points Ai ^ at the moments it £ 
in the same system of coordinants OXYZ, and . 
also the angles of semi-solution .... 

of the diagrams of di recti onali ti es |3 £ , i = 0, /2 


2 . Basic Correlations 

The quantities will be considered, when it is convenient, 
as constitutes of the vector t o All incidental vectors and 

/*£ will be considered uncorrelated between themselves, and their 
universal means as equal to zero: 

£-r 4 f)= 0 , £(A_ri=0, l-57rr:i (') 


where "T" denotes transportation, and "E" dentoes the operation of 
the universal mean. 

We will transfer the origin of the system OXYZ t-o point -^0 and 
designable the system' The radi i -vectors of poi nt: yVf^ 

in the system will equal 


2 



^ = i , 'H-j , 


= r 


O » 


( 2 ) 


and for the covariance matrixes^ of their error^— jn accordance 
with (I), one may write 

(3) 


PkA 




(4) 


It should be noted that the errors^p^, as opposed to the errors 
^ , are correlated,, since ‘ 

-^p/) 

The single vector of direction towards an infinite remote 
source will be exactly the same in both systems, and the radi us -vector 
of a source which is located in a. finite distance will equal 

We will suppose that the velocity of radiation propogation C is 
constant and the wave of radiation for rinfinitely remote sources is 
flat, and for a source located in a finite distance it is spherical. 

Then, in the case of an infinitely remote source the quantities 
T - bound by the correlations 


I 


R^Cy+CX = 0 , 

_i^;i = 1 , 

l 5 co^P; , (■ = r , n- i , 


1 _- , 


( 9 ) 

( 6 ) 


( 7 ) 


where R is the matrix whose columns are = I, n " I, 


In the case of a source within a finite distance, the quantities 
». /pY^ill be bound by the correlations 

( 8 ) 

I V ' “ ' ^ = 1 , Ti - 1 , 

• II 


I Ip* - Pyl - i Pvl = ’ '^ = 1 , 41 - 1 , 

py > py COS , t=0,n-l. (9) 
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Algorithms of the definition and •, with the help of 
systf^ms ,I and II, repsecti vely , will be considered for the cases 
of n = 2, 3, 4, each of these cases we will suppose that 


rank (D) =. r? - < 

where 

D = . 


( 10 ) 

( 11 ) 


If the conditions (10) is not f ul f i 1 1 ed , then it is necessary 
to exclude one. of the ■ poi nts ^ and consider the problem with 
smaller numbers of reference points. 


It should be noted that when n > Z solutions of systems 
I and II can depend on the choice of the original reference point 
M of M* • We will consider, for example, the case where -n = 3. 

. 0 • . t 

Let ij be the rate! of time lag of the transmission of a flare in 
point’j^£ , as compared to the point^-^y , which is determined as a 
result of evaluating the spectra ”and • Depending on the 

choice of th^ initial .point three sets of rates of time lags may be 
obtained: and *, and it and • . 

Since for any method of evaluating spectrT the' condi t'iont^^y = » 

should be fulfilled, only three quantities of T ij can be independent. 
For exampler T^g^nd tg^ , which for convenience we will in- 
dicate as -^(1) , 1 ^^) r a'nd , respectively (if this is com- 
pared to quantities introduced earlier, thent^ = Tg = - T 

If the condition - ... 

a = t + 0 , (^3) 

is fulfilled, then, as is easily seen, the solutions of systems I 
and II do not depend on the choice of the original point. The con- 
dition (13) will be fulfilled, evidently, in the case that the spec- 
tra<SQ S^and Sg are identical. Usualily condition (13) is not ful- 
filled and then the difference of thequantity 0 , from zero will 
characterize the accuracy of the definition of -j^(^)» Another way of 
using condition (13) is by bringing it into the i nterpreta ti on of 
spectra with the definition of this case all spectra 


W1 


as 


11 be processed together. Specifically, ifT^’^'^is defined 
the points of the extremum of several functions (p^(t) > i.e. if 


dz 


= 0 


\z^ 


d'^CDtCZ) , ^ ’ 


It =t^ 


then the condition (13) may be considered by introducing amendments 
into the quanti ti es r ^ It may be shown that the quanti ties , 
which guarantee the attainment of the conditional extremum of func- 
tions ( i . e . by calculating condition (13) ), in a linear approach 

will equal 


+ + ^ = 1 , 2,3 . 

' ’ ' t. 

If al 1 6^ are equal amongst themselves, then 

= , ^= 1 . 2/5 . 


(k) 


(15) 


(16) 


These equations may be applied if 5^ are not know, 


The Case of Two. 


The solution to system I will be the vector , which equals 


^.y ' (9 1 /Pi ) cosa +(u / n) sin a , 


(17) 


where 


coscx , 0 < ot < Jl 


(18) 


and is the random vector, orthogonal Pi ^ 


'H ^ = 0 . 


(19) 


In the celestial sphere on assemble of solutions of form the 

circumference of radius ct. with the center at point, • 
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Errors of measurements A 1 and A , lead to the error ^ ^ of 
the definition of iy which is directed along the normal to this 


circumference at point t. 


The root-mean-square declination CT£ 


of the error A ^depends upon the assumed quantity ^ and equals 

[ 4 ] : ' ^ 


= <^/CPi sin.o, ) , 


( 20 ) 


where 


0-, = (c"x^ + e.^^K J ) 


T 1/2 




( 21 ) 


Since the dimension of the matrix. in this case equals I, is 

simply the dispersion of error • 

The assemble of solutions of system II forms a surface which 
is the hyperboloid of the rotation whose focus lies in points Mq 
and M. The semi-major axis of the rotation equal s |CTj| /2 . When 
; > 0 the solutions are the points of that cavity of the hyper- 

boloid which encompasses the focusMj^and whenT^ 0 the solutions are 
the points of that cavity which encompass the focus M.q* If^ i = 0» 
then the hyperboloid degenerates into a plane which is normal to p j 
and transient through its center. In consequence of the indeter- 
minacy of the solution of system II, it is more expedient to solve 
the problem, which is reversed with respect to the problem of locali- 
zation: the verification of the correspondence to system II of- the cO' 

ordinants of some cosmological objects which could be flare sources. 

If o is the radi us -vector of such a "suspicious" object, then 
a rate of time lag ^; > which corresponds to this object, can be de- 
termined: 




= (lP^-pvj-IPyl)/C.* = I 


( 22 ) 


The hypothesis that a region of space is a "suspected" flare 
source can be assumed by fulfilling the condition 

|r, -i,l /or^ ^ }<s , 
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( 23 ) 


where depends upon the assumed probability of 5 to repudiate 
the correct hypothesis (the level of significance). If the law of 
distribution of the error A is normal, then the quantity /Cj. 
will be determined by the equation 

= 1 -6 . (2 


4. The Case of Three 


In correspondence [4jwith the solution of system I, vector 
will equal 


Cy = {[V. X Pj] ctr^ + ^ Tt ^ c T 2 ^ TL i - tt cp ^ j ^ (25) 


where 


X = D, X 


Pi ^ P 2 y 


(26) 


COS Cp ^ (cx^ D'^ cx.y^^ , 0 < q: ^ Ji/2 . 


(27) 


The expression (27) defines the size of the angle qp between the 
directions towards the source and the reference plane ..which passes 
through points Mq > The condition |coscj/)^i defines 

the region of possible vaUles of T . The region is limited by the 
ellipse 


— r 

cr 


D'^cx - 1 


(28) 


which is written into the square {CTjj ^ and 1^X21^92 • _ On the 

border of the region (28) there exists a single solution iy , which 
lies in the reference plane. Within this area there exist two so- 
lutions of J , which are symmetrical relative to the reference 
plane. For uncovering the non-simplicity it is necessary to appro- 
priate, in this case, condition (7). 
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If ty does not lie in the reference plane, i.e. 

^ 0 , (29) 

then the covariance matrix the error ^^is determined in the 

following manner [^J : 


J(t=S{c^K^*K)S\ (30) 

where S is the matrix of the dimension 3X2, whose columns are 
equal to: 

52= , (31) 

and K is the matrix of dimension 2X2, whose elements equal: 

( 33 ) 

where the covariance ma tri xes X are determined in accordance 
wi th (3) and (4) . 

At least one natural value of the matrix K equal to zero. 
The corresponding natural vector is col linear to ^ ^ . Thus,A"^def i nes 
the dispersion of the plane which is tangent to the celestial 

sphere at point Sy • If in the reference plane, that is 

n = 0 , ( 33 ) 

then the product of the error ^ , which is orthogonal to the re- 

ference plane, turns out to be indeterminant. The root-mean-square 
declination of the projections of error^ ( on the reference plane 
can be determined in the same manner as in the case where /? = 2 

( 20 ). 


For a source located within a finite distance, system II de- 
termines an assemble of possible solutions which form a line of 
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of intersection of two hyperboloids of the rotation. As in the 
case where/? = this instance it is most expedient to solve 

the reversed p^roblem, determining in accordance with (22) the 
quantities Tj and for the coordinants of the supposed 

sources of flares and checking condition (23) for each quantity 

5. The Case of Four 


When measurment errors are lacking the solution of system II 
would be the vector 


o=-RcV. - (34) 

The influence of the measurement errors leads to the failure of 
condition (6). Therefore, in the capacity of £ the normalized 
vector can be used: ^ 



(35) 


However, it is more convenient to use the more general algorithm, 
which determines the solution of system II, and in a particular 
instance the solution of system I also. In accordance with the 
a 1 gori thm L^] the distance to a source 9y is determined from the 
quadratic equation 


In + 5 I = 1 , 


(36) 


P'v = i/pv , 

... 

and the constituents vector equal 

=('P*-c^-D/2 , k 


= i.2,3. 


(37) 

(38) 

(39) 
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Equation (36) has a single positive sol uti on > 0 when 6' < X» i-e., 
if vector'cT i es within the ellipsoid 


cr cr = i . 


(40) 


When S=I, i.e. on the verification of the ellispoid, there exists 
two solutions. One of them corresponds to an infinitely remote source 
( = 0 ). and the other corresponds to a source located at a dis- 

tance 

Py =-P^/2P^S . _ (41) 

t 

When S=I both solutions (36) may be positive. Each solution of H- Y 
determines the direction towards a source at poi nty*^ and equals 


The accuracy of determining the distance to a source depends 
on the position of the source relative to the reference points. 

For the case where the distance to a source significantly exceeds 
the distance between reference points, a simple expression for the 
root-mean-square decl i nati on' CT^ of the error of determining the dis- 
tance can. b e obtai ned : T 

\ ^ ^5 P V / 9 mi-n j 

, . Z n a .1/2 

<^5 - 2 ( 2 r „„) 1 ^ 

where Pmi« , is the least of the distance 9l ’ P 2 » Ps ; ^Tfr,<s.:r i s 

maximum natural value of the matrix;-^^ ; and i s the greatest 

of the maximum natural calues of the matrixes K u . 

i 


If it is demanded that be N times-less than the determined 

distance 9y , then an expression for the greatest distance to 
sources which can be determined with the fullfillment of the follow- 
ing condition can be obtained: 


Pv 


max 


('napceK) 4. 

[ pSrsa.c. ] 




(44) 
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where -A is the rel a ti onshi (T^/’p^^^^iexpres sed i n angular seconds. 
6. Conclusion 


For the localization of infinitely remote sources it is suf- 
ficient to register a flare at three spacecrafts. 

The bivalency arising in this case can be eliminated with the use of 
dilrected detectors. However it is impossible to distinguish a flat- 
wave from a spherical. Therefore, the analysis of registrations of 
flares from nearby sources will determine the direction, which can 
differ considerably from the direction towards the true source. Thus, 
when using three sc it is necessary to verify (the reverse pro- 
blem) the possibility of flare radiation from "suspicious" areas of 
space. The identification of wave form and the subsequent determina- 
tion of the direction towards an infinitely remotesource or coordi- 
nants of nearby sources are possible only with the use of four 
sc. In all cases of localization there occurs the degeneration of 
the measurement system when the source and any two reference points 
are located on one straight line. During the working of the problem 
it is necessary to exclude one of these points and conduct the lo- 
calization using the remaining points. 
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